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ABSTRACT 

In order to test the O star wind-shock scenario for X-ray production in less lu- 
minous stars with weaker winds, we made a pointed 74 ks observation of the nearby 
early B giant, (3 Cru (BO. 5 III), with the Chandra High Energy Transmission Grating 
Spectrometer. We find that the X-ray spectrum is quite soft, with a dominant ther- 
mal component near 3 million K, and that the emission lines are resolved but quite 
narrow, with half- widths of 150 km s . The forbidden-to-intercombination line ratios 
of Ne IX and Mg xi indicate that the hot plasma is distributed in the wind, rather 
than confined near the photosphere. It is difficult to understand the X-ray data in the 
context of the standard wind-shock paradigm for OB stars, primarily because of the 
narrow lines, but also because of the high X-ray production efficiency. A scenario in 
which the bulk of the outer wind is shock heated is broadly consistent with the data, 
but not very well motivated theoretically. It is possible that magnetic channeling could 
explain the X-ray properties, although no field has been detected on (3 Cru. We de- 
tected periodic variability in the hard (hv > 1 keV) X-rays, modulated on the known 
optical period of 4.58 hours, which is the period of the primary (3 Cephei pulsation 
mode for this star. We also have detected, for the first time, an apparent companion 
to (3 Cru at a projected separation of 4". This companion was likely never seen in 
optical images because of the presumed very high contrast between it and (3 Cru in 
the optical. However, the brightness contrast in the X-ray is only 3:1, which is con- 
sistent with the companion being an X-ray active low-mass pre-main-sequence star. 
The companion's X-ray spectrum is relatively hard and variable, as would be expected 
from a post T Tauri star. The age of the (3 Cru system (between 8 and 10 Myr) is 
consistent with this interpretation which, if correct, would add [3 Cru to the roster of 
Lindroos binaries - B stars with low-mass pre-main-sequence companions. 

Key words: stars: early-type - stars: mass-loss - stars: oscillations - stars: pre-main- 
sequence - stars: winds, outflows - stars: individual: f3 Cru - X-rays: stars 



1 INTRODUCTION 

X-ray emission in normal O and early B stars is generally 
thought to arise in shocked regions of hot plasma embed- 
ded in the fast, radiation-driven stellar winds of these very 
luminous objects. High-resolution X-ray spectroscopy of a 
handf ul of bright O star s has basically confirmed this sce- 
nario (|Kahn et all l200ll : iKramer etaU 120031 ; ICohen et ail 
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120061 : iLeutenegger et al.ll2006l h The key diagnostics are line 
profiles, which are Doppler broadened by the wind outflow, 
and the forbidden-to-intercombination (f/i) line ratios in 
helium-like ions of Mg, Si, and S, which are sensitive to the 
distance of the shocked wind plasma from the photosphere. 
Additionally, these O stars typically show rather soft spec- 
tra and little X-ray variability, both of which are consistent 
with the theoretical predictions of the line-driving instabil- 
ity (LDI) scenario l|Owocki et al.l ll988) and in contrast to 
observations of magnetically active coronal sources and mag- 
netically channeled wind sources. 
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Table 1. Properties of /3 Cru 



MK Spectral Type 


BO. 5 III 


Hiltner ct al. (1969'). 


Distance (pc) 


108 ± 7 


Hivvarcos (Perrvman et al. 1997). 


Age (Myrj 


O tO 11 


ijiriaioos (lyooj, cic ocus ct ai. iiyoyj. ocitcin ct ai. (iyy4j ana mis worK. 


6» LD (mas) 


0.722 ±0.023 


Hanbury Brown ct al. (1974). 


M (M ) 


16 


Aerts et al. ( 19981. 


T cS (K) 


27, 000 ± 1000 


Aerts et al. (1998). 


log g (cm s -2 ) 


3.6 ±0.1; 3.8 ±0.1 


Aerts et al. (1998): Mass from Aerts et al. (1998) combined with radius. 


vs'mi (km s — 


35 


Uesuei & Fukuda (1982). 


L (L ) 


3.4 x 10 4 


Aerts et al. (1998). 


R (R ) 


8.4 ±0.6 


From distance (Perrvman et al. 1997) and 9j n Qianbur^Brown ct al. 1974) 


Pulsation periods (hr) 


4.588, 4.028, 4.386, 6.805, 8.618 


Aerts et al. (1998): CuvDers et al. (2002). 


M (M Q yr" 1 ) 


lO- 8 


Theoretical calculation, using Abbott (1982). 


M<? (M yr" 1 ) 


io- n 


Product of mass-loss rate and ionization fraction of C IV (Prinia 1989). 


■Uoo (km s *) 


2000 


Theoretical calculation, using Abbott (198^. 


tJoo (km s -1 ) 


420 


Based on Civ absorption feature blue edge (Prinia 1989). 



In order to investigate the applicability of the LDI wind- 
shock scenario to early B stars, we have obtained a pointed 
Chandra grating observation of a normal early B star, j3 
Crucis (B0. 5 III). The star has a radiation-driven wind but 
one that is much weaker than those of the O stars that 
have been observed with Chandra and XMM-Newton. Its 
mass-loss rate is at least two and more likely three orders 
of magnitude lower than that of the O supergiant £ Pup, 
for example. Our goal is to explore how wind-shock X-ray 
emission changes as one looks to stars with weaker winds, in 
hopes of shedding more light on the wind-shock mechanism 
itself and also on the properties of early B star winds in 
general. We will apply the same diagnostics that have been 
used to analyze the X-ray spectra of O stars: line widths 
and profiles, f/i ratios, temperature and abundance analysis 
from global spectral modeling, and time- variability analysis. 

There are two other properties of this particular star 
that make this observation especially interesting: /3 Cru is 
a P Cephei variable, which will allow us to explore the po- 
tential connection among pulsation, winds, and X-ray emis- 
sion; and the fact that there are several low-mass pre-main- 
sequence stars in the vicinity of (3 Cru. This second fact 
became potentially quite relevant when we unexpectedly dis- 
covered a previously unknown X-ray source four arc seconds 
from P Cru in the Chandra data. 

In the next section, we discuss the properties and en- 
vironment of j3 Cru. In §3 we present the Chandra data. In 
§4 we analyze the spectral and time-variability properties of 
P Cru. In §5 we perform similar analyses of the newly discov- 
ered companion. We discuss the implications of the analyses 
of both P Cru and the companion as well as summarize our 
conclusions in §6. 



a spectral type of B0. 5, this makes P Cru - also known as 
Becrux, Mimosa, and HD 111123 - extremely bright. In fact, 
it is the 19th brightest star in the sky in the V band, and 
as a prominent member of the Southern Cross (Crux), it 
appears in the flags of at least five nations in the Southern 
Hemisphere, including Australia, New Zealand, and Brazil. 

The fundamental properties of P Cru have been exten- 
sively studied and, in our opinion, quite well determined 
now, although binarity and pulsation issues make some of 
this rather tricky. In addition to the Hipparcos distance, 
there is an optical interferometric measurement of the star's 
angular diameter IjHanburv Brown et al.lll974l ). and thus a 
good determi nation of its radius. It is a s pectroscopic binar y 
l|Heintz|[l957l ) with a period of 5 years (|Aerts et al.lll998l V 
Careful analysis of the spectroscopic orbit, in conjunction 
with the brightnes s contrast at 4430 A from interferomet- 
ric measurements (|Popperlll968l ) and comparison to model 
atmospheres, enable determinations of the mass, effectiv e 
temperature, and luminosity of P Cru (|Aerts et al.lll998l ). 
The log g value and mass imply a radius that is consistent 
with that derived from the angular diameter and parallax 
distance, giving us additional confidence that the basic stel- 
lar parameters are now quite well constrained. The projected 
rotational velocity is low (v sin i = 35 km s -1 ), bu t the incli- 
nation is also likely to be low. lAerts et al.l (1 19981 ) constrain 
the orbital inclination to be between 15° and 20°, and bi- 
nary orbital and equ atorial pla nes are well aligned for bina- 
ries with a < 10 AU jHaldl 19941 ). This implies « rot w 120 km 
s -1 and combined with R* = 8.4 Rq, the rotational period 
is Prot ~ 3"?6. The star's properties are summarized in Table 

m 

With Tes = 27, 000 ± 1000 K, P Cru lies near the 
hot edge of the P Cephei pulsation strip in the HR 



2 THE P Cru SYSTEM 

One of the nearest early- type stars to Earth, P Cru is located 
at a distance of only 108 pc l|Perrvman et aLlll997M T l. With 

1 A r ecent re-reduction of the Hipparcos data llvan Leeuwenl 
l2007al fbl) results in a smaller distance of 85_ 6 pc; or 
88^?, pc with the application of the Lutz-Kelker correction 



llMaiz Apellaniz et al.ll2008l) . This re-evaluation of the distance 
has a small effect on the results presented here. It makes the ra- 
dius inferred from the angular diameter correspondingly smaller 
and the gravity larger, as well as decreasing the X-ray luminosity 
inferred from the measured flux. In any case, parallax distance 
determinations may be affected by the spectroscopic binary com- 
panion, and a highly precise distance determination should at- 
tempt to correct for the small orbital angular displacement of 
P Cru during the course of its 5 year orbit. 
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diagram l|Stankov fc Handled 120051 ). Three different non- 
radial pulsation modes with closely spaced periods between 
4.03 and 4.59 hou rs have been identified spectroscopically 
|Aerts et al.l 1 19981 ). Modest radial velocity variations are 
seen at the level of a few km s~ overall, with some indi- 
vidual p ulsation componen ts having somewhat larger am- 
plitudes (|Aerts et al.l If 998h . The three modes are identi- 
fied w ith azimuthal wavenumbers I = 1, 3, and 4, respec- 
tively l|Aerts et al.llf998l : ICuvpers et aLlbOOSft . The observed 
photometric variability is also modest, with amplitudes of 
a few hundredths of a magnitude s een only in the pri- 
mary pulsation mode (|Cuyperslll983l ). More recent, inten- 
sive photometric monitoring with the star tracker aboard 
the WIRE satellite has identified the three spectroscopic 
periods and found two additional very low amplitude (less 
than a millimagnitude) , so mewhat lower frequency compo- 
nents (jCuypers et al.ll2002l ). 

The wind properties of P Cru are not very well known, 
as is often the case for non-sup ergiant B stars. T he stan- 
dard theory of line-driven winds (|Castor et al.lll975l ) (CAK) 
enables one to predict the mass-loss rate and terminal veloc- 
ity of a wind given a line list and stellar properties. There 
has not been much recent work on applying CAK theory 
to B stars, as the problem is actually significantly harder 
than for O stars, because B star winds have many fewer 
constraints from data and the ionization/excitation condi- 
tions are more difficult to calculate accurately. With those 
caveats, howe ver, we have used the CAK parameters from 
lAbbottl (f 982) and calculated the mass-loss rate and termi- 
nal velocity using the st ellar parameters i n Tab le [T] and the 
formalism described by iKudritzki et all (|f 9891 ). The mass- 
loss rate is predicted to be about 10~ 8 Mq yr _1 and the 
terminal velocity, v^, « 2000 km s _1 . There is some uncer- 
tainty based on the assumed ionization balance of helium, 
as well as the uncertainty in the stellar properties (gravity, 
effective temperature, luminosity), and whatever systematic 
errors exist in the line list. The wind models of IVink et al.l 
(2000) also predict a mass-loss rate close to 10 -8 Mq yr _1 
for /3 Cru. 

Interestingly, IUE observations show very weak wind 
signatures - in Si IV and C IV, with no wind signature at all 
seen in Sim. The Civ doublet near 1550 A has a blue edge 
velocity of only 420 km s - , while the blue edge velocity of 
the Si IV line is slightly smaller. The products of the mass- 
loss rate and the ionization fraction of the relevant ion are 
rough ly 10 -11 Mq yr -1 for both of these species (|Prinial 
If 9891 ). Either the wind is much weaker than CAK theory 
(using Abbott's parameters) predicts or the wind has a very 
unusual ionization structure. We will come back to this point 
in §6\ 

P Cru is most likely a member of the Lower Centaurus 
Crux (LCC) subgroup of the Sco-Cen OB associat ion. Al- 
though a classical member of LCC (e.g, Blaauw 1 946f). /3 Cru 
was n ot selected as an LCC member by Ide Zeeuw et ail 
(1999) in their analysis of Hipparcos data. However, the 
timespan of the Hipparcos proper motion measurement (4 
years) is similar to the 5-yr orbital period of the spectro- 
scopic binary companion. The rate of orbital angular dis- 
placement of p Cru is similar to its proper motion, which 
likely influenced the Hipparcos p r oper-m otion measurement. 
For this reason, Ide Zeeuw et al.l (| 1 9991 ) suspect that (3 Cru 
is an LCC member even though it is not formally selected 



as one by their criteria, and indeed iHoogerwerJ (|2000h find 
that f3 Cru is an LCC member based on longer-time-baseline 
proper motion measurements. 

The age of (3 Cru is also roughly consistent with LCC 
membership, given the uncertainty about the age(s) of the 
LCC itself. The age of LCC has been va riously found to 
be from 10-12 Myr Jde Geus et al.l If 9891 ) to 17-23 Myr 
|Mamaiek, Meyer, fc Lieberj 20021 ). and more recent work 
has found evidence for an age spread and spatial s ubstruc- 
ture within the region l|Preibisch fc Mamaiekll2007l ). (3 Cru 
has been claim ed to have an evolutionary age of 8 Myr 
l|Lindroosl If 985) . though the distance was much more un- 
certain when this work was done. Using our adopted stellar 
parameters (Table [T|, f3 Cru has an age of 8-10 Myr whe n 
plotted on the evolutionary tracks of iBertelli et al.l (1 19941 ). 
As such, j3 Cru has evolved off the main sequence, but not 
very far. Its spectroscopic binar y companion is a B2 star 
that is still on the main sequence (|Aerts et al.lll99sl ). and the 
projected X-ray companion discussed below, if it is a phys- 
ical companion, is likely to be a pre-main-sequence (PMS) 
star. There are two other purported wide visual compan- 
ions (P Cru B and C) listed in the Washington Double Star 
catalog (|Worlev fc Douglass! Il997l ). with separations of 44" 
and 370", although t hey are almost certainly not physical 
companions of (3 Cru (|Lindrooslll985l ). 

A small group of X-ray bright pre-main-sequence 
stars within about a degree of P Cru was discovered 
by ROSAT (|Park fc Finlevl Il996l ). These stars are likely 
part of the huge low-ma ss PMS population of the LCC 
|Feigelson fc La wson 1997). They have recently been studied 
in depth (|Alcala et al.ll2002l ). and four of the six have been 
shown to be likely post-T-Tauri stars (and one of those is a 
triple system). They have Ha in emission, strong Li 6708 A 
absorption, and radial velocities consistent with LCC mem- 
bership. Assuming that they are at the same distance as 
P Cru, their ages from PMS evolutionary tracks are 4-10 
Myr. This is somewhat younger than LCC ages derived from 
the more massive members, but this is consistent with a 
larger pattern of systematic age differences for LCC mem- 
bers of different masses an d may indicate systematic e rrors 
in the evolutionary tracks l|Preibisch fc M ama ick [20071 ). 



3 THE Chandra DATA 

We obtained the data we report on in this paper in a single 
long pointing on 28 May 2002, using the Chandra Advanced 
CCD Imaging Spectrometer in spectroscopy mode with the 
High Energy Transmissio n Grating Spectrome ter (ACIS- 
S/HETGS configuration) |Canizares et al .1120051 ). The ACIS 
detector images the dispersed spectra from two grating com- 
ponents of the HETGS, the medium energy grating (MEG) 
and the high energy grating (HEG). The mirror, grating, 
and detector combination has significant response on a wave- 
length range extending roughly from 2 A to 40 A, though 
with a total (including both negative and positive orders) 
effective area of only about 10 cm 2 . The spectral resolution 
exceeds A/AA = 1000 at the long- wavelength end of each 
grating's spectral range, and decreases toward shorter wave- 
lengths. The spatial resolution of Chandra is excellent, with 
the core of the point spread function having a full width of 
less than half an arcsecond. Properties and characteristics of 
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Figure 1. Image of the center of the ACIS detector, showing the 
zeroth order spectra of /3 Cru and its newly discovered companion. 
North is up and east is left. f3 Cru is the source to the northwest. 
Its position in the ACIS detector is coincident with the optical 
position of (3 Cru to within the positional accuracy of Chandra. 
The companion is to the southeast, at a separation of 4.0" and 
with a position angle relative to /3 Cru of 120°. The extraction 
regions for the MEG and HEG (both positive and negative orders) 
are indicated for both sources. It can be seen in this image that the 
position angle is nearly perpendicular to the dispersion directions 
(for both the MEG and HEG) and thus that the dispersed spectra 
of the two sources can be separated. In the on-line color version 
of this figure, hard counts (hv > 1 keV) are blue, soft counts 
(hu < 0.5 keV) are red, and counts with intermediate energies 
are green. It can easily be seen in the color image that the X-ray 
emission of the companion is dramatically harder than that of 
j3 Cru. 



the telescope and d etector are available via the Proposers' 
Observatory Guide (ICXCll2007l ). 

We reran the pipeline reduction tasks in the the 
standard Chandra Interactive Analysis of Observations 
(CIAO) software, version 3.3, using the calibration database, 
CALDB 3.2, and extracted the zeroth order spectrum and 
the dispersed first order MEG and HEG spectra of both 
(3 Cru and the newly discovered companion. The zeroth or- 
der spectrum of a source in the Chandra ACIS-S/HETGS 
is an image created by source photons that pass directly 
through the transmission gratings without being diffracted. 
The ACIS CCD detectors themselves have some inherent 
energy discrimination so a low-resolution spectrum can be 
produced from these zeroth-order counts. 

The zeroth order spectrum of each source has very lit- 
tle pileup (~ 7% for [3 Cru), and the grating spectra have 
none. The HEG spectra and the higher order MEG spectra 
have almost no source counts, so all of the analysis of dis- 
persed spectra we report on here involves just the MEG ±1 
orders. There were typically five to ten pixels across each 
emission line, and few lines have even 100 counts per line. 
The random errors on the source count rates are dominated 
by Poisson photon-counting noise. Thus even the strongest 
lines have signal-to-noise ratios per pixel of only several. The 
background (instrumental and astrophysical) is negligible. 
We made observation-specific response matrix files (rmfs; 
encoding the energy-dependent spectral response) and grat- 
ing ancillary response files (garfs; containing information 
about the wavelength-dependent effective area) and per- 



formed our spectral analysis in both XSPEC 11.3.1 and ISIS 
1.4.8. We performed time- variability analysis on both binned 
light curves and unbinned photon event tables using custom- 
written codes. 

The observation revealed only a small handful of other 
possible sources in the field, which we identified using the 
CIAO task celldetect with the threshold parameter set to 
a value of 3. We extracted background-subtracted source 
counts from three weak sources identified with this proce- 
dure. Their signal-to-noise ratios are each only slightly above 
3. The only other bright source in the field of view is one 
just 4'.'0 to the southeast of (3 Cru at a position angle of 120 
degrees, which is clearly separated from (3 Cru, as can be 
seen in Fig. [1] The position angle uncertainty is less than a 
degree based on the formal errors in the centroids of the two 
sources, and the uncertainty in the separation is less than 
f/.'l. We present the X-ray properties of this source in Sj5]and 
argue in ^that it is most likely a low-mass PMS star in orbit 
around [3 Cru. The B star that dominates the optical light 
of the f3 Cru system and which is the main subject of this 
paper should formally be referred to as f3 Cru A (with spec- 
troscopic binary components 1 and 2). The companion at a 
projected separation of 4" will formally be known as (3 Cru 
D, since designations B and C have already been used for 
the purported companions listed in the Washington Double 
Star Catalog, mentioned in the previous section. Through- 
out this paper, we will refer to (3 Cru A and D simply as 
"/3 Cru" and "the companion." 

We summarize the properties of all the point sources 
on the ACIS chips in Tab. We note that the other three 
sources have X-ray fluxes about two orders of magnitude 
below those of f3 Cru and the newly discovered companion. 
They show no X-ray time variability, and their spectra have 
mean energies of roughly 2 keV and are fit by absorbed 
power laws. These properties, along with the fact that they 
have no counterparts in the 2MASS point source catalog, 
indicate that they are likely AGN. We did not detect any of 
the pre- main-sequence stars seen in the ROSAT observation 
because none of them fell on the ACIS CCDs. We also note 
that there is no X-ray source detected at the location of 
the wide visual companion located 44" from (3 Cru, f3 Cru 
C, which is no longer thou ght to be physically associated 
with f3 Cru (Lindroos 1985). We place a la upper limit of 7 
source counts in an extraction region at the position of this 
companion. This limit corresponds to an X-ray flux more 
than two orders of magnitude below that of j3 Cru or its 
X-ray bright companion. 

The position angle of the newly discovered companion 
is oriented almost exactly perpendicular to the dispersion 
direction of the MEG for this particular observation. This 
can be seen in Fig. [1] where parts of the rectangular ex- 
traction regions for the dispersed spectra are indicated. We 
were thus able to cleanly extract not just the zeroth order 
counts for both sources but also the MEG and HEG spec- 
tra for both sources. We visually inspected histograms of 
count rates versus pixel in the cross-dispersion direction at 
the locations of several lines and verified that there was no 
contamination of one source's spectrum by the other's. We 
analyze these first-order spectra as well as the zeroth-order 
spectra and the timing information for both sources in S|4] 
and g5] 
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Table 2. Point sources in the field 



Name(s) RA^ Dec. total counts_ X-ray flux^, 

(J2000) (J2000) (10~ 14 ergs s-~cra~ 2 ) 



P Cru A; Mimosa; Becrux; HD 111123 


12 


47 


43.35 


-59 


11 


19.2 


3803 ± 12 


186 ± 


1 


P Cru D; CXOU J124743.8-594121 


12 


47 


43.80 


-59 


41 


21.3 


1228 ± 12 


57.3 ± 


.6 


CXOU J124752.5-594345 


12 


47 


52.53 


-59 


1.3 


45.8 


25.9 ±8.6 


1.58 ± 


52 


CXOU J124823.9-593611 


12 


18 


23.91 


-59 


36 


11.8 


47.4 ± 10.4 


3.32 ± 


73 


CXOU J124833.2-593736 


12 


48 


33.22 


-59 


37 


36.8 


29.6 ±9.8 


2.23 ± 


74 



a The positions listed here are not corrected for any errors in the aim point of the telescope. There is an offset of about 0'.'7 
between the position of P Cru in the Chandra observation and its known optical position. This likely represents a systematic 
offset for all the source positions listed here. The statistical uncertainties in the source positions, based on random errors 
in the centroiding of their images in the ACIS detector, are in all cases <0'.'l. 

" Zeroth order spectrum, background subtracted using several regions near the source to sample the background. 
c On the range 0.5 keV < hu < 8 keV and, for the weak sources, based on fits to the zeroth order spectra using a power law 
model with interstellar absorption. For the brighter sources (/3 Cru and its newly discovered companion) the flux is based 
on two-temperature APEC model fits discussed in §4.1 and §5.1. 



4 ANALYSIS OF /3 Cru 
4.1 Spectral Analysis 

The dispersed spectrum of j3 Cru is very soft and relatively 
weak. We show the co-added negative and positive first- 
order MEG spectra in Fig. [2] with the strong lines labeled. 
The softness is apparent, for example, in the relative weak- 
ness of the Nex Lyman-a line at 12.13 A compared to the 
Neix He a complex near 13.5 A, and indeed the lack of any 
strong lines at wavelengths shorter than 12 A. 



4-1.1 Global thermal spectral modeling 

We characterized the temperature and abundance distribu- 
tions in the hot plasma on j3 Cru by fitting Astrophysi- 
cal Plasma Emission Code (APEC) thermal, o ptically thin, 
equil ibrium plasma spectral emission models (jSmith et al.l 
l200ll ) to both the zeroth and first order spectra. We included 
a model for pileup in the fit to the zeroth order spectrum 
l|Davisll200ll ). For the dispersed MEG spectrum, we adap- 
tively smoothed the data for the fitting in ISIS. This effec- 
tively combines bins with very few counts to improve the 
statistics in the continuum. We also modified the vAPEC 
("v" for variable abundances) implementation in ISIS to al- 
low the helium-like forbidden-to-intercombination line ratios 
to be free parameters (the so-called Q parameter, which is 
the ratio of the forbidden plus intercombination line fluxes 
to the flux in the resonance line, Q = was also allowed 
to be a free parameter). We discuss these line complexes 
and their diagnostic power later in this subsection. The in- 
clusion of these altered line ratios in the global modeling 
was intended at this stage simply to prevent these particular 
lines from influencing the fits in ways that are not physically 
meaningful. 

The number of free parameters in a global thermal 
model like APEC can quickly proliferate. We wanted to 
allow for a non-isothermal plasma, non-solar abundances, 
and thermal and turbulent line broadening (in addition to 
the non-standard forbidden-to-intercombination line ratios) . 
However, we were careful to introduce new model param- 
eters only when they were justified. So, for example, we 
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Figure 2. The MEG spectrum of P Cru, with negative and posi- 
tive first orders coadded and emission lines identified by ion. Bins 
are 5 mA wide, and each emission line spans several bins. We 
show a smoothed version of these same data in Fig. [4] 



only allowed non-solar abundances for species that have sev- 
eral emission lines in the MEG spectrum. And we found we 
could not rely solely on the automated fitting procedures in 
XSPEC and ISIS because information about specific physi- 
cal parameters is sometimes contained in only a small por- 
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tion of the spectrum, so that it has negligible effect on the 
global fit statistic. 

We thus followed an iterative procedure in which we 
first fit the zeroth-order spectrum with a solar-abundance, 
two-temperature APEC model. We were able to achieve a 
good fit with this simple type of model, but when we com- 
pared the best fit model to the dispersed spectrum, system- 
atic deviations between the model and data were apparent. 
Simply varying all the free model parameters to mi nimize 
the fit statistic (we used both the Cash C statistic ( Cash ! 
1 19791 ) and x 2 1 periodically comparing the results for the two 
statistics) was not productive because, for example, the con- 
tinuum in the MEG spectrum contains many more bins than 
do the lines, but the lines generally contain more indepen- 
dent information. And for some quantities, like individual 
abundances, only a small portion of the spectrum shows 
any dependence on that parameter, so we would hold other 
model parameters fixed and fit specific parameters on a re- 
stricted subset of the data. Once a fit was achieved, we would 
hold that parameter fixed and free others and refit the data 
as a whole. 

Ultimately, we were able to find a model that provided 
a good fit to the combined data, both zeroth-order and dis- 
persed spectra. The model parameters are listed in Tab. [3] 
and the model is plotted along with the zeroth order spec- 
trum in Fig.[3]and the adaptively smoothed first order MEG 
spectrum in Fig. [3] There are some systematic deviations 
between the model and the data, but overall the fit qual- 
ity is good, with xl = 1-09 for the combined data. The 
formal confidence limits on the fit parameters are unrealis- 
tically narrow, at least partly because the statistical error 
cannot be fully accounted for. For example, the abundance 
determinations are each made on narrow wavelength inter- 
vals and thus the temperature cannot reliably be allowed 
to vary along with the abundances in the confidence limit 
testing procedure. We return to this point below when we 
discuss the abundance determinations. The overall result, 
that a thermal model with a dominant temperature near 3 
million K and an emission measure of ~ 3.5 x 10 53 cm -3 and 
abundances that overall are slightly sub-solar is required to 
fit the data, is quite robust. 

The model we present provides a better fit for the first- 
order MEG spectrum, which has 443 bins when we adap- 
tively smooth it. The model gives xt = 1-00 for the first 
order MEG spectrum alone, while it gives xi = 1-42 for the 
zeroth order spectrum alone, which has 90 bins. We suspect 
the relatively poor agreement between the CCD spectrum 
and the MEG spectrum stems from the small, but non-zero, 
transmission of red optical light through the MEG and HEG 
grati ngs and through the ACIS optical blocking filter l|Wolkl 
I2002T ). Being so bright, (3 Cru is more susceptible to this ef- 
fect than nearly any other object. We examined the Chandra 
ACIS spectrum of Vega, which is very bright in the optical 
and X-ray dark. Its CCD spectrum, which is presumably 
due entirely to the optical contamination, looks very similar 
to the excess above the model seen in the f3 Cru spectrum. 
The optical contamination is particularly noticeable in the 
0.22 — 0.5 keV ACIS-S band, which we have chosen not to 
model. There appears to be some contamination up to ~ 0.7 
keV, however. 

In Fig. [4] where we show our global best-fit model su- 
perimposed on the adaptively smoothed MEG spectrum, the 
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Figure 3. The best-fit two-temperature, variable abundance 
thermal spectral (APEC) model (dotted line), corrected for 
pilcup, superimposed on the zeroth-order spectrum of f3 Cru. The 
model shown is the best-fit model fit to both the zeroth and first 
order spectra, with parameters listed in Tab. \3\ It seems likely 
that some of the systematic deviations in the softer channels are 
due to optical contamination, as discussed in §4.1.1. 



continuum is seen to be well reproduced and the strengths of 
all but the very weakest lines are reproduced to well within a 
factor of two. It is interesting that no more than two temper- 
ature components are needed to fit the data; the continuum 
shape and line strengths are simultaneously explained by 
a strong 2.5 million K component and a much weaker 6.5 
million K component. 

The constraints on abundances from X-ray emission line 
spectra are entangled with temperature effects. So, it is im- 
portant to have line strength measurements from multiple 
ionization stages. This is strictly only true in the MEG spec- 
trum for O and Ne, but important lines of Mg, Si, and Fe 
from non-dominant ionization stages would appear in the 
spectrum if they were strong, so the non-detections of lines 
of Mg XII, Si xiv, and Fe xviii provide information about 
the abundances of these elements. These five elements are 
the only ones for which we report abundances derived from 
the model fitting (see Tab.|3|. We allowed the nitrogen abun- 
dance to be a free parameter of the fit as well, but do not 
report its value, as it is constrained by only one detected line 
in the data. The uncertainties on these abundances are dif- 
ficult to reliably quantify (because of the interplay with the 
temperature distribution) but they are probably less than 
a factor of two, based on our experience varying these pa- 
rameter values by hand. They are certainly somewhat larger 
than the formal confidence limits listed in Tab. [3] but even 
given that, there is an overall trend of slightly subsolar abun- 
dances, consistent with a value perhaps somewhat above half 
solar. Only the oxygen abundance seems to be robustly and 
significantly sub-solar. We note, however, that in general, 
the abundances w e find are significantly closer to solar than 
those reported bv lZhekov fe Pallal (|2007l ). 

In all fits we included ISM absorption assuming a 
hydrogen column density of 3.5 x 10 19 cm' 2 . This value 
is co nsistent with the very low E(B — V) = 0.002 ± 
.011 (Nic mczura fe Daszvnska-Daszkiewiczl 120051 ) and is 
the value determine d for the nearby star HD 110956 
l|Fruscione et al . 1994). This low ISM column produces very 



Chandra spectroscopy of (3 Cru 7 




Wavelength (A) 



Ne IX 



Ne IX 



NeX 




Fe XVII 



Ne IX 



FeXVII - 




13 

Wavelength (A) 



VIII 



FeXVII 




Ju 




18 



Wavelength (A) 



VII 



J 



N VII 



22 23 
Wavelength (A) 



Figure 4. The best-fit two-temperature, variable abundance thermal spectral (vAPEC) model (red, dotted line; parameters listed in Tab. 
\3\ same model as shown in Fig. [3}, superimposed on the adaptively smoothed co-added MEG first-order spectrum of f3 Cru. Emission lines 
are labeled and color coded according to element, with the abundances of each element listed at the top of the figure. Some non-detected 
lines are labeled, including the forbidden line of O VII near 22.1 A and the Lyman-alpha line of Si XIV near 6.2 A. 



little attenuation, with only a few percent of the flux ab- 
sorbed even at the longest Chandra wavelengths. The over- 
all X-ray flux of j3 Cru between 0.5 keV and 8 keV is 
1.86 ± .01 x 10~ 12 ergs s" 1 cm" 2 , implying L x = 2.6 x 10 30 
ergs s _1 , which corresponds to logL x /I/Boi = —7.72. 



4-1.2 Individual emission lines: strengths and widths 

Even given the extreme softness of the dispersed spec- 
trum shown in Figs. [2] [3] and [4] the most striking thing 
about these data is the narrowness of the X-ray emission 
lines. Broad emission lines in grating spectra of hot stars 
are the hallmark of wind-shock X-ray emission dKahn et al 



200ll : ICassin"elh" et al 1 l200ll : iKramer et al]|2003l : ICohen et al 
2006). While the stellar wind of f3 Cru is weaker than that 
of O stars, it is expected that the same wind-shock mech- 



anism that operates in O stars is also responsible for the 
X-ray emission in these slightly later-type early B stars 
jCohen et alJll997h . The rather narrow X-ray emission lines 
we see in the MEG spectrum of Cru would thus seem to 
pose a severe challenge to the application of the wind-shock 
scenario to this early B star. We will quantify these line 
widths here, and discuss their implications in Sj6] 

We first fit each emission line in the MEG spectrum in- 
dividually, fitting the negative and positive first order spec- 
tra simultaneously (but not coadded) with a Gaussian pro- 
file model plus a power law to represent the continuum. In 
general, we fit the continuum near a line separately to estab- 
lish the continuum level, using a power-law index of n — 2, 
(F v oc v n , which gives a flat spectrum in F\ vs. A). We then 
fixed the continuum level and fit the Gaussian model to the 
region of the spectrum containing the line. We used the C 
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Table 4. Emission Lines in the HETGS/MEG Spectrum 



ion 


A lab 


\ b 

A obs _ 


shift 


half-width^ 


normalization 


corrected flux " 






(A) 


(mA) 


(mA) " 


(10 -6 ph s" 1 cm" 2 ) 


(10 — 15 erg s — 1 cm -2 ) 


Si xm 


6.6479 (r) 






1 —1.7 


1.06+°« 
— 0.34 


3.15+:' 2 ? 

— 1.01 


Si xm 


6.6882 (i) 






1.7 


3l + !!:2 

J — 0.17 


93+?,-!? 

u,i,J - 0.51 


Si xm 


6.7403 (f) 






1.7 


77 +0 ,'Si 3 . 

u 1 1 -0.26 


2 26+°'?,I 
^ -0.76 


Me xi 


9.1687 (r) 






1 5+?' 2 

± -° — 1.5 


^'^ -0.52 


4.81+HI 

1.13 


Me xi 


9.2312 (i) 






1.5 


-, 9r -+0.50 
-'-■'''-'-0.40 


2 70+i-°! 

' "—0.86 


Me xi 


9.3143 (f) 






1.5 


82+°- 4 ? 


1 56+i' ,? 

1,J "-0 .71 


Ne IX 


11.5440 






1 3+J 2 , 8 

1.3 


6 0+H 


10 5 +2 .'? 

lu,J -2.5 


Ne X 


12.1339 


12.1331 ± 0.0019 


-0.8 ± 1.9 


6.3+ 2 i 

— 2.3 


13.0 +2 g 


21.5+ 4 '? 

— 3.1 


Ne IX 


13.4473 (r) 






7.0 ± 1.3 


39.9+ 4 '° 

—5.2 


59.6+«'° 

— 6.6 


Ne IX 


13.5531 (i) 






7.0 


37 4+ 4 - 2 

Jl • -4.9 




Ne IX 


13.6990 (f) 






7.0 


3.4+H 

— 1.5 


4.9+?? 

— 2.2 


Fs xvn 


15.0140 


15 0146+ 0021 
lO.U14D_ 0020 


6+ 21 

U.O_ 2 g 


10 0+ 2 ' 6 
1U.U_ 2 


0U -°-ll.l 


81 1+ 7 ' 6 


Fe xvn 


15.2610 






i -i 7+3.7 
' -3.0 


28.1+f; 2 


37.1+s-; 


O viii 


16.0059 


16.0053+°;°^ 




9+ 51 


32.2+^ 




Fe xvm 


16.0710 








S-2+25 


10.3+11 


Fe xvn 


16.7800 


16 7756+ 0044 
1D '" OD -0.0010 


_ 4 4+4-4 




31.8+*-? 


38.3_ g ' 9 


Fe xvn 


17.0510 


17.04771°;°°?° 


-3.31 2 ;? 


s.8±i:? 


67.3+lj 


78.4+9'° 


Fe xvn 


17.0960 


i™938lS:S2S 


z,z -1.6 


8.8 


53-2±|;i 


62.9+ 7 ' 4 


O vii 


18.6270 






3 9+ 12 ' 4 


21-115.9 


22 9+ 9 '° 
ZZ - 9 -6.4 


O viii 


18.9689 


18.967X1°;°°" 




10.4 ± 1.3 


260l 2 4 


278+J 2 


O vii 


21.6015 (r) 


21 6016+ 0018 

Z1.DU1D_ Q017 




11.0 ± 1.8 


295 ± 34 


279 ± 30 


O vii 


21.8038 (i) 


91 7QS9+0.0012 

zj.. iaoz_ OQOg 




■= 4 +l.6 


414l 42 


389+J9 


O vi£ 


22.0977 (f) 








< 19.3 


< 17.9 


N vii 


24.7810 






7 7+3.7 
'■'-3.2 


92 ±21 


75+ 18 
'°-16 



a Lab wavelength taken from the Astrophysical Plasma Emission Database (APED) llSmith et al. 2001). For closely spaced doublets 
(e.g. Lyman-alpha lines) an emissivity-weighted mean wavelength is reported. The resonance (r), intercombination (i), and forbidden 
(f) lines of Si, Mg, Ne, and O are indicated. 

^ If no value is given, then this parameter was held constant at the laboratory value when the fit was performed. 

c For closely spaced lines, such as the He-like complexes, we tied the line widths together but treated the common line width as a free 
parameter. This is indicated by identical half-width values and formal uncertainties listed on just the line with the shortest wavelength. 
** Corrected for an ISM column density of Ng = 3.5 X 10 19 cm~ 2 using Morrison & McCammon (1983) photoelectric cross sections 
as implemented in the XSPEC model wabs. 

e The normalization and corrected flux values for the forbidden line of O VII are based on the 68% confidence limit for a Gaussian line 
profile with a fixed centroid and a fixed width (of 7 mA). The data are consistent with a non-detection of this line. 



Table 3. Best-fit parameters for the two-temperature APEC 
thermal equilibrium spectral model 



parameter 


value 


90% confidence limits 


Ti (K) 


2.5 x 10 6 


2.4 : 2.6 x 10 6 


SMi (cm" 3 ) 


3.3 x 10 53 


3.1 : 3.5 X 10 53 


T 2 (K) 


6.5 x 10 6 


6.2 : 6.8 X 10 6 


£M 2 (cm" 3 ) 


2.8 x 10 52 


2.3 : 3.3 X 10 52 


Dturb (km s _1 ) 


180 


90:191 


o 


0.35 


0.30:0.40 


Ne 


0.78 


0.70:0.99 


Mg 


0.52 


0.29:0.62 


Si 


1.05 


0.68:1.56 


Fe 


0.80 


0.70:0.90 



The abunda nces are the (linear) fractio n of the solar abundance, 
according to lAnders fe Grevessd l ll989h . 



of individual lines. We report the results of these fits, and 
the derived properties of the lines, in Tab. [4] and show the 
results for the line widths in Fig. [5] Note that the half- width 
at half-maximum line width, for the higher signal-to-noise 
lines (darker symbols in Fig. [5}, is about 200 km s~ . The 
weighted mean HWHM of all the lines is 143 km s" 1 , which 
is indicated in the figure by the dashed line. The typical ther- 
mal width for these lines is expected to be roughly 50 km 
s _1 , so thermal broadening cannot account for the observed 
line widths. 

To explicitly test for the presence of broadening beyond 
thermal broadening, when we fit the two-temperature vari- 
able abundance APEC models to the entire MEG spectrum 
we accounted explicitly for turbulent and thermal broad- 
ening. We found a turbulent velocity of 180 km s - , corre- 
sponding to a half- width at half-maximum turbulent veloc- 
ity of 150 km s _1 , from a single global fit to the spectrum. 



statistic (|Cashlll979h . which is appropriate for data where 
at least some bins have very few counts, for all the fitting 



2 As implemented in ISIS, the Gaussian line profile half- width is 
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We will characterize the overall line widths as fhwhn 



150 



km s~ for the rest of the paper. The thermal velocity for 
each line follows from the fitted temperature(s) and the 
atomic mass of the parent ion. This turbulent broadening 
is consistent with what we find from the individual line fit- 
ting, as shown in Fig. [S] We also find a shift of —20 km s _1 
from a single global fit to the spectrum. We note that this 
20 km s _1 blue shift of the line centroids is not significant, 
given the uncer tainty in the absolute wavelength calibration 
of the HETGS ^Marshall et alj|2004h . 

Because the emission lines are only barely resolved, 
there is not a lot of additional information that can be ex- 
tracted directly from their profile shapes. However, because 
the working assumption is that this line emission arises in 
the stellar wind, it makes sense to fit line profile models 
that are specific to stellar wind X-ray emission. We thus fit 
the empi rical wind-profile X-ray emission line model devel- 
oped by lOwocki fc Cohenl (2001) to each line, again with 
a power-law continuum model included. This wind-profile 
model, though informed by numerical simulations of line- 
driving instability (LDI) wind shocks ( Owocki et al.l 1 19881 ; 
ICohen et al.lll99 6; Feld meier et al.lll997 ). is phenomenolog- 
ical, and only depends on the spatial distribution of X-ray 
emitting plasma, its assumed kinematics (described by a 
beta velocity law in a spherically symmetric wind), and the 
optical depth of the bulk cool wind. For the low-density wind 
of (3 Cru, it is safe to assume that the wind is optically thin 
to X-rays. And once a velocity law is specified, the only free 
parameter of this wind-profile model - aside from the nor- 
malization - is the minimum radius, -R m i n , below which there 
is assumed to be no X-ray emission. Above i? m i n , the X-ray 
emission is assumed to scale with the square of the mean 
wind density. 

When we fit these wind-profile models, we had to de- 
cide what velocity law to use for the model. From CAK 
ijCastor et al.l Il975h theory, we expect a terminal velocity 
of roughly 2000 km s _1 . The f3 parameter of the standard 
wind velocity law governs the wind acceleration according 
to v = Voo(l — R./r)' 3 so that large values of f3 give more 
gradual accelerations. Typically (3 ~ 1. When we fit the 
wind-profile model with Vaa = 2000 km s" 1 and (3 — 1 to 
the stronger lines, we could only fit the data if i? m i n < 1-1 
R» . This puts nearly all of the emission at the slow-moving 
base of the wind (because of the density-squared weighting 
of X-ray emissivity and the v _1 r~ 2 dependence of density), 
which is unrealistic for any type of wind-shock scenario, but 
which is the only way to produce the relatively narrow pro- 
files in a model with a large wind terminal velocity. This 
quantitative result verifies the qualitative impression that 
the relatively narrow emission lines in the MEG spectrum 
are inconsistent with a standard wind-shock scenario in the 
context of a fast wind. 

However, although the theoretical expectation is for a 
fast wind with a terminal velocity of roughly 2000 km s , 
the UV wind line profiles observed with IUE tell a different 
story. The UV line with the strongest wind signature is the 
C IV doublet at 1548,155 1 A, which sh ows a blue edge veloc- 
ity of only 420 km s" 1 (|Prinialll989r ). We show this line in 
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Figure 5. Half-width at half-maximum line widths derived from 
fitting Gaussian line-profile models to each individual line in the 
MEG first-order spectrum. The strongest, well behaved, isolated 
lines are indicated with darker symbols. The lower signal-to- 
noise, blended, or otherwise less reliable lines are shown in gray. 
68% confidence limit error bars are shown for all the lines. The 
weighted mean HWHM velocity for all the lines is indicated by 
the dashed horizontal line, while the HWHM that corresponds 
to the turbulent velocity derived from the global spectral fit is 
indicated by the dotted horizontal line (recall that v tuI ^, = 180 
km s _1 corresponds to a half width at half maximum of 150 km 
B- 1 ). 
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given by 0.83 yVyj + v tmb' so HWHM that corresponds to 
the best-fit turbulent velocity value is 150 km s _1 . 
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Figure 6. The Civ resonance doublet from several coadded 
archived IUE observations. A velocity scale (km s — 1 ) in the frame 
of the star is shown above the data. These lines show the char- 
acteristic blue-shifted absorption that is expected from a stellar 
wind, but the lines are not only weak compared to theoretical 
expectations but also much less wind broadened. 



Fig. [6] The Si ill line shows no wind broadening and the Si 
IV line shows a slightly weaker wind signature than the C IV 
feature. Certainly the terminal velocity of the wind could be 
significantly higher than this, with the outer wind density 
being too low to cause noticeable absorption at high veloci- 
ties. However, because the empirical evidence from these UV 
wind features indicates that the terminal velocity of the wind 
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Table 5. Helium-like complexes 



km 


n = f/i± 


Auv (At 


H v (ergs s 1 cm 2 Hz 


Rfir (R*ll 


Rmin (R*l^ 


c - (/+») 


Ne ix 
Mg xi 


0.09 ± .05 
0.65 ± .44 


1248, 1273 
998, 1034 


1.35 x 10~ 3 
1.39 x 10~ 3 


2 9+ - 7 
, 9 +1.4 


l.37±;" 
1.80t£; 33 


1.05 
1.01 



a From Gaussian fitting. The values we obtain from the APEC global fitting within ISIS are consistent with the f/i values reported here. 
^ The UV wavelengths at which photoexcitation out of the upper level of the forbidden line occurs. 

c The assumed photospheric Eddington flux at the relevant UV wavelengths, based on the TLUSTY model atmosphere, which we use to 
calculate the dependence of f/i on the radius of formation, -Rfl r . 

Formation radius using eqn. (lc) of lBlumenthal et al.l jl972ft and eqn. (3) of iLeutenegger et al. 

I (I2006T) . and based on the f/i values in 

the second column. 

e From wind profile fitting, assuming (3 = 1 and letting Voa be a free parameter. The Q values in the last column also come from the profile 
fitting. 
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Figure 7. The Fe XVII line at 15.014 A with three different wind 
profile models. The red, dashed line is a constant-velocity wind 
model with Voo = 280 km s — 1 (the best-fit value for a constant 
outflow velocity, which produces emission lines with Uhwhm ~ 
150 km s" 1 ). The model with Vac = 2000 km s _1 , /3 = 1, and 
^min = 1-5 R* is represented by the blue, dot-dash line. An 
infinitely narrow model is shown as the green, dotted line. The 
residuals for each model fit are shown in the lower panel, as red 
circles for the global best-fit, modestly broadened (voa = 280 km 
s _1 ) model, green squares for the narrow profile model, and blue 
triangles for the broad wind model. The wind model with the 
higher velocity (v^c = 2000 km s _1 ) clearly does not provide 
a good fit, while the narrower constant velocity (voc = 280 km 
s -1 ) wind model does. And while the very narrow model cannot 
be absolutely ruled out, the Voo = 280 km s _1 model is preferred 
over it with a high degree of significance. All the models shown 
here have been convolved with the instrumental response function 
(rmf) and multiplied by the wavelength-dependent effective area 
(90-rf). 



may actually be quite low, we refit the strong lines in the 
MEG spectrum with the wind-profile model, but this time 
using a terminal velocity of 420 km s _1 . The fits were again 
statistically good, and for these models, the fitted values of 
Rmka were generally between 1.3 R* and 1.5 R„, which are 
much more reasonable values for the onset radius of wind- 
shock emission. 

We also considered a constant- velocity model (/3 = 0), 
where the X-rays perhaps are produced in some sort of ter- 
mination shock. For these fits, J? m in hardly affects the line 



width, so we fixed it at 1.5 R* and let the (constant) terminal 
velocity, Voo, be a free parameter. For the seven lines or line 
complexes with high enough signal-to-noise to make this fit- 
ting meaningful, we found best-fit velocities between 270 km 
s _1 and 300 km s _1 for five of them, with one having a best- 
fit value lower than this range and one line having a best-fit 
value higher than this range. The average wind velocity for 
these seven lines is = 280 km s _1 . All in all, the emission 
line widths are consistent^ with a constant-velocity outflow 
of 280 km s _1 , and are inconsistent with a wind-shock origin 
in a wind with a terminal velocity of 2000 km s" 1 . The lines 
are narrower than would be expected even for a wind with a 
terminal velocity not much greater than the observed C IV 
blue edge velocity, v = 420 km s _1 (based on our /3 = 1 fits 
as well as the f3 = fits). The X-ray emitting plasma must 
be moving significantly slower than the bulk wind. In Fig. [7] 
we show the best-fit constant velocity model (woo = 280 km 
s _1 ) superimposed on the Fe xvn line at 15.014 A, along 
with the higher velocity model based on the theoretically ex- 
pected terminal velocity, which clearly does not fit the data 
well. We also compare the best-fit wind-profile model to a 
completely narrow profile in this figure. The modestly wind- 
broadened profile that provides the best fit is preferred over 
the narrow profile at the 99% confidence level, based on the 
C stati stic values o f the r espective fits (Numerical Recipes, 
Ch. 14; |Press et all (|l982T l). 

4-1.3 Helium-like forbidden-to-intercombination line 
ratios 

The final spectral diagnostic we employ is the UV-sensitive 
forbidden-to-intercombination line ratio of helium-like ions. 
These two transitions to the ground-state have vastly dif- 
ferent lifetimes, so photoexcitation out of the ls2s 3 Si 
state, which is the upper level of the forbidden line, to the 

3 The measured line widths (of Uhwhm ~ 150 km s — 1 derived 
from both the Gaussian fitting and the global thermal spectral 
modeling) are expected to be somewhat smaller than the modeled 
wind velocity, as some of the wind, in a spherically symmetric 
outflow, will be moving tangentially to the observer's line of sight. 
So, a spherically symmetric, constant velocity wind with v^, = 
280 km s _1 is completely consistent with emission lines having 
HWHM values of 150 km s" 1 . 
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Figure 8. The Ne IX resonance-intercombination-forbidden com- 
plex in the co-added MEG data, with the best-fit three-Gaussian 
model overplotted (top). The constraints on the measured f/i 
ratio (0.09 ± .05), from the fit shown in the top panel, are in- 
dicated by the cross-hatched area's intersection with the y-axis 
(bottom). The solid black curve is the model for the line ratio, as 
a function of radius. The cross-hatched region's intersection with 
the model defines the 68% confidence limit on r/R* of (1.9:3.6). 
The dashed line represents the best-fit value, f/i = 0.09, and the 
corresponding i?fl r = 2.9 R». 



Is2p Pi, 2 state, which is the upper level of the intercombi- 
nation line, can decrease the forbidden-to-intercombination 
line ratio, 1Z = f/i. Because the photoexcitation rate de- 
pends on the local UV mean intensity, the f/i ratio is a di- 
agnostic of the distance of the X-ray emitting plasma from 
the photosphere. This diagnostic has been applied to many 
of the hot stars that have Chandra grating spectra, and typ- 
ically shows a source location of a t least half a stellar radius 
above the photosphere. Recently, iLeutenegger et al.l (2006) 
have shown that the f/i ratios from a spatially distributed 
X-ray emitting plasma can be accounted for self-consistently 
with the line profile shapes in four O stars, providing a uni- 
fied picture of wind-shock X-ray emission on O stars. 

The lower-Z species like oxygen do not provide any 
meaningful constraints, as is the case for O stars, and there 
are not enough counts in the Si xiii complex to put any 
meaningful constraints on the source location. Only the Ne 
IX and Mg XI complexes provide useful constraints. For both 
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Figure 9. The Mg XI resonance-intercombination-forbidden com- 
plex in the co-added MEG data, with the best-fit wind-profile 
models overplotted (top). The constraints on the formation ra- 
dius, i?fi r , are shown in the bottom panel. 



Ne IX and Mg XI we find Rr t w 3± 1R» . Of course, the line- 
emitting plasma is actually distr ibuted in the wind, and 
so we also fit the same modified lOwocki fc Cohenl (|200ll ) 
wind profile model that ILeutenegger et al.l ( 20061 ) used to 
fit a distributed X-ray source model to the He-like com- 
plexes, accounting simultaneously for line ratio variations 
and wind-broadened line profiles. With the wind velocity 
law and Q = £^ as free parameters, strong constraints 
could not be put on the models. Typically, these distributed 
models of f/i gave onset radii, Rmi n , that were quite a bit 
lower than Ra T . This is not surprising, because a distributed 
model will inevitably have some emission closer to the star 
than i?fl r to compensate for some emission arising at larger 
radii. We list the results of such a wind-profile fit to the 
He- like complexes in Tab. [5] along with the R^ v values. The 
two helium-like complexes and the best-fit models are shown 
in Figs. [5] and [5] along with the modeling of the f/i ratios 
and the characteristic formation radii, Ra T . We note that 
the data rule out plasma right at the photosphere only at 
the 68% confidence level. 

We derive a single characteristic radius of formation, 
i?fi r , for each complex by calculating the dependence of 
the f/i ratio on the distance from the photosphere, using 
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Figure 10. Rejection probability (significance level increasing 
downward) of assumed period, according to the Kuiper statis- 
tic. The five known optical pulsation periods llAerts et al.lll99sl ; 
ICuvpers et~aT]|2002fi are indicated by the vertical lines. The pri- 
mary pulsation period is the third of these. 

a TLUSTY model atm osphere (lLanz fc Hubeny||2007h . We 
used the formalism of iBlumenthal et al.l (I1972T ). equation 
(lc), ignoring collisional excitation out of the 3 Si level, and 
explicitly expressing the dependence on the dilution factor, 
and th us the radius, as in equation (3) in lLeutenegger et al.l 
( 2006). The measurement of f/i from the data can then be 
used to constrain the characteristic formation radius, Ra T . 

4.2 Time Variability Analysis 

We first examined the Chandra data for overall, stochastic 
variability. Both a K-S test applied to the unbinned photon 
arrival times and a x 2 fitting of a constant source model 
to the binned light curve showed no evidence for variability 
in the combined zeroth order and first order spectra. This 
is to be expected from wind-shock X-rays, which generally 
show no significant variability. This is usually interpreted as 
evidence that the X-ray emitting plasma is dis tributed over 
nume rous spatial regions in the stellar wind (|Cohen et al.l 
1 19971 ). Separate tests of the hard (hv > 1 keV) and soft 
(hv < 1 keV) counts showed no evidence for variability in 
the soft data, but some evidence for variability (K-S statistic 
implies 98% significance) in the hard data. 

Because /3 Cru is a /3 Cephei variable, we also tested the 
data for periodic variability. Both the changes in the star's 
effective temperature with phase and the possible effects of 
wave leakage of the pulsations into the stellar wind could 
cause the X-ray emission, if it arises in wind shocks, to show 
a dep endence on the pulsation phase (|Owocki fc Cranmerl 
2002). Our Chandra observation covers about four pulsation 
periods. To test for periodic varia bility, we used a variant 
of the K-S test, the Kuiper test (|Paltanil [20o3 ). which is 
also applied to unbinned photon arrival times. To test the 
significance of a given periodicity, we converted the photon 
arrival times to phase and formed a cumulative distribution 
of arrival phases from the events table. We then calculated 
the Kuiper statistic and its significance for that particular 
period. By repeating this process for a grid of test periods, 
we identified significant periodicities in the data. 

When we applied this procedure to all the zeroth-order 
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Figure 11. The light curve of the hard counts in the zeroth 
and first order spectra of the primary, along with the best-fit 
sinusoidal model. 

and first-order counts we did not find any significant peri- 
odicities. However, applying it only to the hard counts pro- 
duced a significant peak near the primary and tertiary op- 
tical p ulsation periods (f i and fs) of 4.588 hours and 4.386 
hours (|Aerts et alj Il998l ) , as shown in Fig. [10] This peak 
is significant at the 99.95% level. It is relatively broad as 
the data set covers only about four cycles of the pulsation, 
so it is consistent with bo th pulsation mode s, but not with 
the secondary (/i) mode (|Aerts et al.| [l998). No significant 
periods are found in the soft bandpass, despite the higher 
signal-to-noise there. To quantify the level of periodic vari- 
ability we found in the hard counts, we made a binned light 
curve from these counts and fit a sine wave with a period of 
4.588 hours, leaving only the phase and the amplitude as free 
parameters. We find a best-fit amplitude of 18%. We show 
the light curve and this fit in Fig. [11] and note that there 
appears to be some additional variability signal beyond the 
strictly periodic component. We also note that there is a 
phase shift of about a quarter period (A(f> = 0.27 ± .04) be- 
tween the times of maximum X-ray and optical light, with 
the X-rays lagging behind the optical variability. To make 
this assessment, we compared the time of maximum X-ray 
light from the sinusoidal fitting shown in Fig.[TT]to the time 
of maximum optical light in the WIRE data, based on just 
the primary pulsation mode (Cuypers 2007, private commu- 
nication). We note that <j> = 0.25 corresponds to the maxi- 
mum blue shift in the radi al velocity variability of /3 Cephei 
stars l|Mathias et al.lll992T ). 



5 ANALYSIS OF THE COMPANION 

We subjected the newly discovered companion to most of 
the same analyses we have applied to /3 Cru, as described 
in the previous section. The exceptions are the wind-profile 
fitting, which is not relevant for the unresolved lines in the 
companion's spectrum, and the tests for periodic variability, 
since none is expected. 



Our working hypothesis is that the X-ray bright com- 
panion is a low-mass pre-main-sequ ence star, similar to 



those found in the ROSAT pointing i 


Park & Finlevlll996l; 


Fcigclson & Lawson 1997: Alcala ct al. 


2002). In the follow- 
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Figure 13. The best-fit two-temperature, variable abundance thermal spectral (APEC) model (green, dotted line), superimposed on 
the adaptively smoothed co-added MEG first-order spectrum of the companion (black, solid line). 
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Figure 12. The zeroth order spectrum of the companion (open 
squares) is shown along with that of /3 Cru (filled circles). The 
companion's spectrum is significantly harder, and is well-fit by 
a two-temperature thermal model with component temperatures 
several times higher than those found for ft Cru. The best-fit 
models - to the combined zeroth order and grating spectra - arc 
shown as the dashed (/3 Cru) and dotted (companion) lines. 

ing two subsections, we marshal evidence from the spectral 
and time- variability properties to address this hypothesis. 

5.1 Spectral Analysis 

The spectrum of the companion is significantly harder than 
that of P Cru. This can be seen in both the zeroth-order 
spectrum (see Fig. [12] and Fig. [1} and the dispersed spec- 
trum (see Fig. I13p . The dispersed spectrum has quite poor 
signal-to-noise, as the hardness leads to fewer counts per 
unit energy flux. 

We fit a two-temperature optically thin coronal equi- 
librium spectral model to the zeroth order spectrum of 
the companion. As we did for the primary, we compared 
the best-fit model's predictions to the MEG spectrum and 



found good agreement. This source is significantly harder 
than P Cru, as indicated by the model fit, which has com- 
ponent temperatures of 6.2 and 24 million K. The emis- 
sion measures of the two components are 1.8 x 10 52 cm" 3 
and 4.2 x 10 52 cm -3 , respectively, yielding an X-ray flux 
of 5.71 ± .06 x 1CT 13 ergs s" 1 cm" 2 , corrected for ISM ab- 
sorption using the same column density that we used for 
the analysis of the spectrum of P Cru. This flux implies 
L x = 8.0 x 10 29 ergs s -1 . No significant constraints can be 
put on the abundances of the X-ray emitting plasma from fit- 
ting the Chandra data. This best-fit two-temperature APEC 
model is shown superimposed on the adaptively smoothed, 
coadded MEG spectrum in Fig. [13] (it is also shown, as the 
dotted histogram, in Fig. I12p . These MEG data are much 
noisier than those from P Cru, but the overall fit is rea- 
sonable, and the stronger continuum, indicative of higher 
temperature thermal emission is quite well fit. We note that 
only three of the line complexes have even a dozen counts 
in this spectrum. 

Unlike P Cru, the companion has lines that are narrow - 
unresolved at the resolution limit of the HETGS. This, and 
the hardness of the spectrum, indicates that the companion 
is consistent with being an active low-mass main-sequence 
or pre-main-sequence star. 

The f/i ratios that provided information about the UV 
mean intensity and thus the source distance from the pho- 
tosphere in the case of P Cru can be used in coronal sources 
as a density diagnostic. The depopulation of the upper level 
of the forbidden line will be driven by collisions, rather than 
photoexcitation, in the case of cooler stars. Some classical T 
Tauri stars (cTTSs) have shown anomalously low f/i ratios 
compared to mai n sequence stars and e ven weak-line T Tau ri 
stars (wTTSs) l|Kastner et al.l |2002| ; ISchmitt et all 120051 ). 
The best constraint in the companion's spectrum is provided 
by Ne IX, which has f/i = 2.2 ± 1.1. This value is consistent 
with the low-density limit of f/i — 3.0 (|Porquet fc Dubaul 
2000), which is typically seen in wTTSs and active main 
sequence stars but is sometimes altered in cTTSs. Using 



14 D. Cohen et al. 



0.06 
0.05 



<j 0.03 



0.01 
0.00 



1(1 It 



2x10 4x10 6x10 8x10 

Time (s) 

Figure 14. The binned light curve of the companion, including 
all zcroth- and first-order counts. 




Figure 15. The hardness ratio as a function of time, defined 



as HR 



CR(>lkcV)-CR(<lkcV) 



where CR denotes the count 



CR(>lkcV) + Cfl(<lkcV) ' 
rate in the combined zeroth order and first order (MEG) data, 
over the indicated photon energy ranges. We fit a linear function 
to the data, and find x 2 = 25 for 24 degrees of freedom for a 
constant model (slope of zero), so a constant source model cannot 
be rejected outright. However, the best-fit slope of —1.9 X 10~ 6 
cts s — 2 (solid line) gives x 2 = 20, for 23 degrees of freedom, which 
is a statistically significant improvement over the constant model. 



the calculations in iPorauet fc Dubaul (|2000D . we can place 
an upper limit on the electron density of n c « 10 12 cm -3 , 
based on the la lower limit of the f/i ratio. 

5.2 Time Variability Analysis 

The companion is clearly much more variable than is /3 Cru. 
The K-S test of the unbinned photon arrival times compared 
to a constant source model shows significant variability in all 
energy bands, though the significance is higher in the hard 
X-rays (hv > 1 keV). The binned light curve also shows 
highly significant variability, again with the null hypothesis 
of a constant source rejected for the overall data, and both 
the soft and hard counts, separately. An overall decrease in 
the count rate is clearly seen as the observation progressed 

(Fig. EH. 

The short rise followed by a steady decline of nearly 



50% in the X-ray count rate over the remainder of the ob- 
servation is qualitatively consistent with fla re activity seen 
in late-type PMS stars (Fav ata et al.l 12005). although the 
rise is not as rapid as would be expected for a single, large 
flare. We examined the spectral evolution of the companion 
during this decay, and there is mild evidence for a softening 
of the spectrum at late times, but not the dramatic spec- 
tral evolution expected from a cooling magnetic loop after 
a flare. For example, the hardness ratio (the difference be- 
tween the count rate above 1 keV and below 1 keV, normal- 
ized by the total count rate) is shown in Fig. 1151 There is 
a mild preference for a negative slope over a constant hard- 
ness ratio (significant at the 95% level), but with the data 
as noisy as they are, a constant hardness ratio cannot be 
definitively ruled out. 



6 DISCUSSION AND CONCLUSIONS 

Although it is very hot and luminous, (3 Cru has a Chandra 
spectrum that looks qualitatively quite different from those 
of O stars, with their significantly Doppler broadened X-ray 
emission lines. Specifically, the emission lines of f3 Cru are 
quite narrow, although they are resolved in the MEG, with 
half- widths of ~ 150 km s _1 . These relatively narrow lines 
are incompatible with the wind-shock scenario that applies 
to O stars if the wind of j3 Cru has a terminal velocity close 
to the expected value of w 2000 km s _1 . They are, how- 
ever, consistent with a spherically symmetric outflow with 
a constant velocity of Voo = 280 km s _1 . 

If the wind terminal velocity is significantly less than 
2000 km s _1 - perhaps just a few hundred km s _1 above 
the 420 km s _1 maximum absorption velocity shown by the 
strongest of the observed UV wind lines - a wind-shock ori- 
gin of the X-ray emission is possible, but still difficult to 
explain physically. One wind-shock scenario that we have 
considered relies on the fact that the low-density wind of 
P Cru should take a long time to cool once it is shock 
heated. If a large fraction of the accelerating wind flow 
passes through a relatively strong shock front at r ~ 1.5 
R», where the pre-shock wind velocity may be somewhat 
less than 1000 km s _1 , it could be decelerated to just a few 
100 km s _1 , explaining the rather narrow X-ray emission 
lines. Similar onset radii are seen observationally in O stars 
llKramer et all 120031 ; ICohen et all 120061 ; iLeutenegger et al.1 
120061 ). The hottest plasma component from the global ther- 
mal spectral modeling requires a shock jump velocity of 
roughly 700 km s _1 , with the dominant cooler component 
requiring a shock jump velocity of roughly 400 km s _1 . 

The f/i ratios of Ne IX and Mg XI indicate that the 
hot plasma is either several stellar radii from the photo- 
sphere or distributed throughout the wind starting at sev- 
eral tenths of a stellar radius. This is broadly consistent with 
this scenario of an X-ray emitting outer wind. The lack of 
overall X-ray variability would also be broadly consistent 
with this quasi-steady-state wind-shock scenario. And the 
periodic variability seen in the hard X-rays could indicate 
that the shock front conditions are responding to the stellar 
pulsations, while the lack of soft X-ray variability would be 
explained if the softer X-rays come from cooling post-shock 
plasma in which the variability signal has been washed out. 

The overall level of X-ray emission is quite high given 
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the modest wind density, which is a problem that has been 
long-recognized in early B stars l|Cohen et al.lll997h . A large 
fraction of the wind is required to explain the X-ray emis- 
sion. Elaborating on the scenario of a shock-heated and in- 
efficiently cooled outer wind, we can crudely relate the ob- 
served X-ray emission measure to the uncertain mass-loss 
rate. 

Assuming a spherically symmetric wind in which the 
X-rays arise in a constant-velocity post-shock region that 
begins at a location, Ro, somewhat above the photosphere; 
say Ro w 1.5 (in units of R*), we can calculate the emission 
measure available for X-ray production from 

f f°° 2 

SM = I n a n-ndV = / n c nn'lnr dr. 
J J R 

Again assuming spherical symmetry and a smooth flow with 
a constant velocity for r > Ro, we have 



SM = 1.3 x 10" 



Mi 



(cm" 3 ) 



for the emission measure of the entire wind above r = Ro, 
where M_g is the mass-loss rate in units of 10 _9 Mq yr _1 
and vioo is the wind velocity in units of 100 km s . Equat- 
ing this estimate of the available emission measure to that 
which is actually observed in the data, 3.5 x 10 53 cm -3 , we 
get 



M_o = 



SM 



1.3 x 10 54 



—Rovioo 



for Ro w 1.5 and v — 280 km s _1 (vioo = 2.8), which is the 
velocity obtained from fitting the line widths in the MEG 
spectrum. Thus a mass-loss rate of order 10 _9 Mq yx _1 is 
consistent with the observed X-ray emission measurqj. This 
type of simple analysis relating the observed X-ray emis- 
sion measure in early B stars to their mass-loss rate in 



4 We hasten to point out that our analysis ignores the possibility 
that less than the entire wind volume beyond r = Ro is X-ray 
emitting. Certainly a filling factor term could multiply the mass- 
loss rate in the expression for the emission measure. Ignoring such 
a finite filling factor would lead us to underestimate the mass-loss 
rate. However, another simplification we have made here involves 
the assumption of a smooth, spherically symmetric distribution. 
Any clumping, or deviation from smoothness, in the post-shock 
region will lead to more emission measure for a given mass-loss 
rate (and assumed Ro and vioo)- Ignoring this effect will lead us 
to overestimate the mass-loss rate. Thus the mass-loss rate we 
calculate here of M fs 10~ 9 Mq yr _1 should be taken to be a 
crude estimate, subject to a fair amount of uncertainty. However, 
the two major oversimplifications in this calculation will tend 
to cancel each other out, as we have just described. We thus 
conclude that a mass- loss rate of order 10 — 9 Mq yr~ 1 is sufficient 
to explain the observed X-ray emission, assuming a constant, slow 
outflow of the post-shock plasma above an onset radius of about 
1.5 R*. This mass-loss rate is lower than the predictions of CAK 
theory but significantly larger than the observed mass loss in the 
C IV UV resonance line, implying that the actual mass-loss rate of 
f) Cru is in fact lower than CAK theory predicts, but that some of 
the weakness of the observed UV wind lines is due to ionization 
effects. Of course, the presence of a large X-ray emitting volume 
in the outer wind provides a ready source of ionizing photons that 
can easily penetrate back into the inner, cool wind and boost the 
ionization state of metals in the wind. 



the context of shock heating and inefficient radiative cool- 
ing in the outer wind was applied to R OSAT observation s 
of the early B giants e and j3 CMa by iDrew et all l|l994h . 
who, like we have here, found lower mass-loss rates than 
expected from theory. A downward revision in the mass- 
loss rates of early B stars might not be surprising given 
recent work indicating t hat the mass-loss rates of O stars 
have been overestimated (jBouret et al ]|2005l ; lFullerton~ et al.l 
l2006t ICohen et alj|2007h . 



Given the estimated mass-loss rate of M 



10" 



Mq yr -1 for (3 Cru, we can assess the cooling time, or cool- 
ing length, of the purported wind shocks. A crude estimate 
of the cooling time and length of shock-heated wind mate- 
rial can be made by comparing the thermal energy content 
of the post-shock plasma, 

3 i 
-nkT » 0.1 (ergs cm ) 



E- 



5 x 10" 



I -1 -3s 

(ergs s cm ) 



to the radiative cooling rate, 
dE 

— — = n 

dt 

These numbers come from taking expected values of the 
wind density based on the assumptions of spherical sym- 
metry of the wind and a constant outflow velocity, 



x 10 5 



M- 



(cm 



R 2 vioo 

where R is the radial location in units of stellar radii. Using 
Vioo = 2.8 and M_g = 1, the density is just n » 10 8 cm -3 at 
R = 1.5. We use temperatures of T = 2.5 and 6.5 x 10 6 K, 
given by the spectral fitting discussed in §4.1, assuming that 
some of the softer thermal component is produced directly 
by shock heating (rather than by radiative cooling of hotter 
plasma) since the emission measure of the hot component 
is so small. We take the integrated line emissivity to be 
A ~ 5 x 10~ 23 ergs s" 1 c m 3 for a plasma of several million 
K (|Beniamin et al.l l200ll l . The characteristic cooling time 
thus derived, 



^cool — 



3kT 
2nA 



1 to 2 x 10° 



(8), 



is many times longer than the characteristic flow time, 
t flow EE - « — » 2 x 10 4 (s). 

V V 

Thus, material shock heated to several million degrees at 
half a stellar radius above the surface will essentially never 
cool back down to the ambient wind temperature by radia- 
tive cooling, motivating the picture we have presented of 
a wind with an inner, cool acceleration zone and a quasi- 
steady-state outer region that is hot, ionized, X-ray emit- 
ting, and moving at a more-or-less constant velocity (be- 
cause it is too highly ionized to be effectively radiatively 
driven) . 

Although the scenario we have outlined above is phe- 
nomenologically plausible, there are several significant prob- 
lems with it. First of all, the behavior of wind shocks pro- 
duced by the line-driving instability is not generally as we 
have described here. Such shocks tend to propagate out- 
war d at an appreciable fraction of the ambie nt wind veloc- 
ity (jOwocki et al.lll988l : lFeldmeier et alJll997h . Whereas the 
scenario we have outlined involves a strong shock that is 
nearly stationary in the frame of the star, as one would get 
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from running the wind flow into a wall. There is, however, 
no wall for the wind to run into. It is conceivable that a 
given shock front forms near r — 1.5 R* and propagates 
outward, with another shock forming behind it (upstream) 
at roughly the same radius, so that the ensemble of shocks 
is close to steady state. However, the narrow observed X- 
ray lines severely limit the velocity of the shock front, as 
the post-shock velocity in the star's frame must be greater 
the shock front's velocity. Yet the observed lines are con- 
sistent with a post-shock flow of only v = 280 km s _1 . A 
related, second, problem with the scenario is that the shocks 
required to heat the plasma to T ~ 2.5 to 6.5 x 10 6 K are rel- 
atively strong for a weak wind. Smaller velocity dispersions 
are seen in statistical analyses of the LDI-indu ced struc- 
ture in O star winds (|Owocki fe Runacresl I2002T ). Perhaps 
the purported wind shocks could be seeded by pulsations 
at the base of the wind. Seeding the instability at the base 
has been sho wn to lead to strong er shocks and more X- 
ray emission l|Feldmeier et al.ll 19971 ). though the shocks still 
propagate relatively rapidly away from the star. 

Another problem with this scenario is the very large X- 
ray production efficiency that is implied. Again, this general 
problem - that early B stars produce a lot of X-rays given 
their weak winds - has been known for quite some time. 
Quantifying it here, using the outer wind shock scenario pre- 
sented above, and the mass-loss rate of M « 10 -9 Mq yr _1 
implied by this analysis, an appreciable fraction of the wind 
material is heated to X-ray emitting temperatures at any 
given time. And since the wind must be nearly completely 
decelerated when it passes through the shock front, a similar 
fraction of the available wind kinetic energy would go into 
heating the wind to T > 10 6 K. Th e self-excited LDI typ i- 
cally is much less efficient than this (jFeldmeier et al.lll997h . 

If this scenario of a nearly completely shocked-heated 
outer wind with a low velocity is not likely, then what are 
the alternatives? As already stated, a more standard model 
of embedded wind shocks formed by the LDI and outflowing 
with the wind cannot explain the very modest X-ray emis- 
sion line widths nor the relatively large emission measures 
(or shock heating efficiency). However, a standard coronal- 
type scenario, as applied successfully to cool stars, cannot 
explain the data either. Coronal sources do not have emis- 
sion lines that can be resolved in the Chandra HETGS. And 
the f/i ratios would be lower than observed if the hot plasma 
were magnetically confined near the surface of the star. And 
of course, there is no reasonable expectation of a strong dy- 
namo in OB stars. The relative softness and lack of variabil- 
ity would also seem to argue against an origin for the X-rays 
in a traditional corona at the base of the wind. 

Perhaps magnetic fields are involved, but just not in the 
context of a dynamo and corona. Magnetic channeling of 
hot-star winds has been shown to efficiently produce X-rays 
and t hat X-ray emission has only modestly broad emission 
lines l|Gagne et alj|2005h . There is, however, no indication 
that /3 Cru has a magnetic field. However, a field of only 
B ~ 5 G would be sufficiently strong to channel the low- 
density wind of f3 Cru. The degree of magnetic confinement 
and channeling can be described by the parameter 

B 2 Rl 

MVoo 

where values of 77, > 1 imply strong magnetic confinement 



and B is the equatorial field stre ngth for an assumed dipole 
field (|ud-Doula fc Owockil l2002h . We note that a stronger 
field has a lready been detecte d in another /3 Cephei star, (3 
Cep itself (|Donati et al.ll2002h . though attemp ts to detect a 
field on /3 Cru have not yet been successful (|Hubrig et alj 
2006). The upper limits on the field strength are roughly 
an order of magnitude higher than what would be needed 
for significant confinement and channeling, though. We also 
note that the radius of th e last closed loop scales as r? / 5 
|ud-Doula fc Owockil liooll ). so a field of 50 G - still unde- 
tectable - would lead to a magnetosphere with closed field 
lines out to 3 R*, consistent with the value found for i?fi r in 
§4.1.3. 

In a more speculative scenario, it is possible that the 
material giving rise to the X-ray emission is near the base 
of the wind, and is falling back toward the star - a failed 
wind. This scenario has been invoked to expla in the X-ray 
emission from r Sco (B0.2 V) (|Howk et alj200d ). which sub- 
se quently had a comp lex surface magnetic field detected on 
it (|Donati et al J 120061 ) . If material gets accelerated off the 
surface by radiation pressure, but after moving a distance 
off the surface finds that it can no longer be accelerated - 
perhaps because of the formation of optically thick clumps - 
those clumps may stall or even fall back toward the star and 
interact with the ambient wind, leading to shock heating of 
wind material and material on the surfaces of the clumps. 
The problems with this scenario include the ad hoc nature of 
the clump formation and wind fall-back (though this behav- 
ior is in fact seen in conjunction with magne tic channeling 
(lud-Doula fc Owock i 2002; Gagne et al. I l2005l )) and the f/i 
ratios which indicate that the hot material is several stel- 
lar radii above the surface. Chandra observations of r Sco 
do show f/i ratios consistent with hot plasma significantly 
above the photosphere and they also show only modestly 
broadened X-ray emission lines, similar to what we see in 
/ 3 Cru, although th e spectrum of r Sco is significantly harder 
|Cohen et al.ll2003l ). 

Finally, perhaps the physical characteristics of a wind 
with a slow, hot outer component that we described in our 
initial scenario is correct - that is what the f/i ratios and 
the line widths are telling us - but the heating mechanism 
is not LDI-related wind shocks. Instead a more complete 
plasma treatment may be required to describe at least the 
outer wind, where the densities are below 10 8 cm -3 . In such 
low density winds, the radiatively driven ions can decouple 
from the protons and rapidly accelerate, leading to frictional 
heating dS pringmann & Pauldrach Il992l; iGavlev fc Owockil 
ll994l ; lKrticka fc Kubatll200ll ; IOwocki fc Pulsll2002l ). Heating 
to X-ray emitting temper atures is predicted to o ccur only 
at about spectral type B5 (|Krticka fc Kubat]|200ll ). but per- 
haps if the wind mass-loss rates of B stars are lower than 
anticipated significant frictional heating could occur in early 
B stars, like [i Cru. 

We point out that the energy requirements for heating 
the low density (n < 10 8 cm -3 ) outer wind of fi Cru are not 
severe at all. The heating mechanism might even involve 
wave propagation if there is even a weak magnetic field on 
this star. And once significantly heated, such low density 
circumstellar matter will tend to remain h ot, due to the 
inefficiency of radiative cooling l|Owockill2004h . In this sense, 
the hot plasma in the wind of /3 Cru could be considered a 
corona, and though likely initiated by radiation pressure as 
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described in CAK theory, the low-density, slow, far wind of 
f3 Cru might even be driven by gas pressure gradients. 

None of these scenarios of X-ray production on j3 Cru 
are both physically natural and also in agreement with all 
the available data. Although, if a magnetic field were to be 
detected on (3 Cru, then the magnetically channeled wind 
shock interpretation would be quite reasonable. However, we 
stress that the main characteristics that constrain any mod- 
els that may be put forward are quite secure: the plasma 
temperature is several million degrees; the emission is not 
highly variable, except for the periodic hard X-ray variabil- 
ity; the X-ray emission lines are broadened but only mod- 
estly so; and the X-ray production efficiency (if the X-ray 
emission is related to the star's wind) is rather high. Finally, 
we note that the periodic signal in the hard X-rays and its 
phasing with the radial velocity variability might be a clue 
to the X-ray production mechanism. 

Put in the context of the empirical X-ray spectral prop- 
erties of OB stars, the case of j3 Cru may be extreme but 
it seems to lie at one end of a continuum of line-profile 
behavior in hot stars with r adiation-driven st ellar winds. 
Th e late- type O s tars 5 Ori llMiller et al.ll2002h and a Ori 
A |Skinner et al.1 12004 120081 ) have lines that are broad, 
but not as broad as would be expected in the standard 
wind-shock scenario given their known wind properties. 
Perhaps there is a gradual transition from stars with X- 
ray emission line s consistent with their fast stellar winds , 
such as C Pup ilKahn et al l l200ll; ICassinelli et al.l l200ll ; 
iKramer et 81.1120031 ) and C Ori l|Cohen et al.ll2006l ). through 
some late O stars with X-ray lines of intermediate width, 
to early B stars like (3 Cru for which Chandra can barely 
resolve their lines. 

The newly discovered companion to /3 Cru is most 
likely a low-mass pre-main-sequence star, sim ilar to the post 
T Tauri stars found in the ROSAT data (|Park fc Finlevl 
1 19961 ; lAlcala et ail [20 02). It has a hard, thermal spectrum 
as PMS stars do, and shows signif icant variability, again, as 
PMS stars do (jFavata et al.ll2005l ). The age of f3 Cru is less 
than the p re-ma in-sequence lifetimes of most low-mass stars. 
iLindroosl (1 19861 ) identified a few dozen B stars with likely 
PMS companions. If this newly discovered companion in- 
deed is a PMS star, then /3 Cru would be a Lindroos binary. 
As such, it would be useful for testing evolutionary models 
if an optical spectrum of the companion can be obtained. 
It would be especially helpful if enough information about 
the mid-B spectroscopic com panion (which is supposed to 
still be on the main sequence (|Aerts et al.lll998l )) can also be 
obtained. There are systematic age differences as a function 
of spectral ty pe in the LCC which seem t o reflect problems 
with models llPreibisch fc Mamaiekll2007h . and the analysis 
of the three coeval stars in the /3 Cru system could provide 
important constraints on the evolutionary models. 

Given the companion's X-ray luminosity of nearly 10 30 
ergs s" 1 cm -2 at the assumed distance of 108 pc, even a 
G star somewhat above the main sequence would have a 
fx/ fv ratio consistent with those seen in PMS stars. The 
nearby K and M post-T Tauri stars detected in the ROSAT 



jointing are typically 



13 and have fx/fv 



10" 



typic ally 

(|Park fc Finlevl Il996l ). If the companion has a similar X- 



ray-to-V-magnitude flux ratio, then it is a 12 th magnitude 
PMS K star. If the magnitude difference between /3 Cru and 
this newly discovered companion is roughly 11, then at a 



separation of 4" and with the primary being so bright, de- 
tecting it with a ground-based optical or IR telescope will be 
a challenge, though such data would obviously be very useful 
for confirming its PMS status and characterizing its proper- 
ties. Any accretion disk it might have would be significantly 
irradiated by UV light from the primary. Finally, given a 
projected separation of 430 AU, the companion should have 
an orbital period (around /3 Cru and its spectroscopic com- 
panion) of at least 1000 years. 
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